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Abstract 

We study Bose-Einstein gas for an arbitrary power low interaction C a r~ a . 
This is done by the Hartree Fock Bogoliubov (HFB) approach at T < T c and 
the mean field approach at T > T c . Especially, we investigate the ground 
state property of Bose gas interacting through the Van der Waals — C$r~ 6 
plus Csr~ 3 interactions. We show that the ground state under this interaction 
is stable if the ratio of coupling constants is larger than that of the critical 
curve. We find that the C^r~ z term plays an important role for the stability 
of the ground state when the density of atoms becomes sufficiently large 
at low temperature. Further, using the numerical values of C3 and Cq, we 
confirm that the ground state of alkali atoms are stable. 
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1 Introduction 



Bose-Einstein condensation (BEC) was first observed in dilute ultracold al- 
kali atoms of rubidium |IJ, lithium [FJ] and sodium ||. Furthermore, Mewes 
et al. and Ensher et al. measured the condensate fraction and the energy 
of rubidium atoms 0, 0. It is shown that the transition temperature T is 
shifted not only by the interaction effect but also by the finite size effect by 
a few percent, where T denotes the transition temperature of the noninter- 
acting Bose gas within the external field in the thermodynamic limit. This 
means that both the interaction and the finite size effects play an important 
role in the real Bose gases. In fact, the number of trapped atoms is typically 
N < 10 7 , and this may not be sufficiently large to take the thermodynamic 
limit. 

On the other hand, these experimental findings have stimulated much interest 
in the theory of the interacting Bose gas. Since BEC occurs when atoms are 
dilute and cold, we can treat the interaction of atoms as the two-body inter- 
action. In this case, the s-wave scattering length characterizes the strength 
of the two-body interaction. Under these conditions, the Gross Pitaevskii 
(GP) equation can well describe the behavior of the interacting Bose gas at 
zero temperature |6], |7j . This is a mean field approach for the order parame- 
ter associated with the condensate. Using the GP equation, several authors 
studied the ground state and the excitation properties of the condensate || . 
To study BEC at finite temperature, the GP equation at zero temperature 
was extended by Griffin M. It is called the Hartree Fock Bogoliubov (HFB) 
theory. In particular, Popov approximation to the HFB theory has been 



employed to explain the experimental results ]T0|. 

In the realistic point of view, BEC experiments are carried out in magnetic 
traps. In this situation, spin-polarized atoms interact through a triplet po- 
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tential. As long as atoms remain polarized, they cannot form molecules. For 
alkali atoms, the triplet potential has many bound states which allow them 
to recombine into molecules. Since this recombination can only occur in a 
three-body scattering, it cannot occur for sufficiently low density of atoms. 
Thus, the two-body scattering is dominant. Therefore, spin-polarized atoms 
can remain the gas through dipole two-body scattering which flip the spin 
to untrapped state, and then atoms can produce BEC. In this sense, it is of 
particular importance to find a way to investigate BEC with more realistic 
interactions. 

In this paper, we study BEC for an arbitrary power low interaction C a r~~ a 
which is more realistic than the ordinary contact interaction (5-function). 
This interaction plays an important role for alkali atoms ||11||. The atomic 



interaction V of alkali atoms can be approximately written as |12 



V = V c + V d + V hf + V z + V so , (1) 

where V c is the central force of the interaction. Further, V c can be written 
in terms of the electron exchange interaction V ex and the dispersion force 
— C$r~ 6 — Cgr -8 — Ciqt^ 10 ■ ■ ■ . The magnetic-dipole interaction Vd and the 
hyperfine interaction Vhf behave like C^r' 3 , where the coupling constant C3 
is the spin part of the interaction. The last two terms, Vz and V so represent 
the Zeeman interaction and the spin-orbit interaction, respectively. These 
interactions are important at high density. 

Here, we study Bose gas interacting through the Van der Waals —C^r^ 6 plus 
C%r~ 3 interactions. We show that the ground state under this interaction is 
stable if the ratio of the coupling constants is larger than that of the critical 
curve. We find that the C%r~ z term plays an important role to stabilize the 
ground state when the density of atoms becomes large at low temperature. 
Using the numerical values of C3 and C$, we confirm that the ground state 



of alkali atoms are stable. 

This paper is organized in the following way. In the next section, we derive 
the effective Hamiltonian for the interacting Bose gas. Then, in section 
3, we study the ground state stability for the Van der Waals plus C^r -3 
interactions. Section 4 summarizes what we have clarified in this paper. 

2 The Effective Hamiltonian for Interacting 
Bose Gas 

In this section, we derive the effective Hamiltonian for the interacting Bose 
gas. Since Bose gas dramatically changes its behavior below the transition 
temperature T c , we derive the effective Hamiltonian with the HFB approach 
at T < T c . On the other hand, we derive the effective Hamiltonian with the 
mean field approach at T > T c . 

The Hamiltonian for the interacting Bose gas confined in a harmonic oscil- 
lator potential can be written as 

h 2 



H = J dr¥{r) 



V 2 H — mujr 2 

2m 2 



+1 I dr 1 dr 2 ^t(r 1 )^(r 2 )V(|r 1 - r 2 |)*(r 2 )*( ri ), (2) 
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where ^(r) is the boson field operator. The two-body atomic interaction 
V(|ri — r 2 |) is given as 

V(\n-r 2 \) = (3) 

Fl — r 2 



Now, we write down the second quantized Hamiltonian for Eq. (2) |13| . The 
boson field operator for the ideal system can be expanded by plane waves, 
but for the general case, the corresponding field operator is a sum over all 
normal modes Jl4] 



$(r) = ^2a uX u(r), (4) 



V 



when x^(r)'s are any complete set of normalized single-particle wave function, 
and a v is a bosonic annihilation operator for the single-particle state v. From 
Eq. (4), the Hamiltonian can be written as 



(5) 



vXv'X' 



with 



V v>\' 



-s 



■h 2 + 



rfr 1 rfr 2 xt(ri)xl(r 2 )l / (|r 1 - r 2 |)xy(r 2 )xy( r i)- 



(6) 
(7) 



Here, we expand V(|ri — r 2 |) in terms of the Legendre polynomial Pi 

oo 

T^dn — r 2 j) = ^u,(ri,r2)P,(cos0i 2 ), 

1=0 

oo I 

= EE 

1=0 m=-l 

where f;(ri,r 2 ) is given by 



47T 



2Z + 1 



vi(r u r 2 )Y l *(e 1 ,ip 1 )Y l , m (e 2 ,i P2 ), (8) 



21 4- 1 Z" 1 

«i(ri,r 2 ) = y dtV(|n-r 2 |)fl(t). 



(9) 



In the case of the power low potential V(|ri — r 2 |) = C a |ri — r 2 | a , we can 
integrate vi(ri, r 2 ); 



^(n,?" 2 ) = 
for ri > r 2 , and 



r(l)r(< + f ) r\ (a a 1 3 4 

r(i + i)r(f)^ 21 U + ' 2 _ 2' ' + o' 3 ) ' ( 10 ) 



2'r? 



T(/ + ±)r(f> 2 ^ 



2 2' 



for r 2 > ri, where 2 -F\(a, 6, c, x) denotes hypergeo metric function. Thus, we 
obtain the second quantized Hamiltonian for Bose gas interacting through 
the power low interaction. From now on, we present the mean field Hamil- 
tonian and the HFB Hamiltonian. 
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2.1 Mean Field Theory 



First, we present the mean field theory for the interacting Bose gas |T3 



This theory is particularly valid for T > T c But for T < T c , this 



theory is inadequate, since the low energy part of the excited states plays an 
important role at low temperature. This effective Hamiltonian H e ff is given 
by the diagonal part of the Hamiltonian (5), and can be written as 



Heff = ^ T vv h v + ^ v vv >h v h v >, (12) 



vv' 



where n u = ala u and 



v w = j 't/r^lx^rOlVdri-rsDIx,^)! 2 , (13) 

for v = v 1 and 

J dv 1 dv 2 \ Xv {v 1 )\ 2 V{\v 1 -v 2 \)\xA^)\ 2 
+ / dr 1 dr 2 xI(r 1 )xU r 2)^ / (|ri - r 2 |)x„(r 2 )x^(ri) (14) 



v uv > 



for v ^ v' . We rewrite Eq. (12) in the following way 

H eff = H + H', (15) 



with 



H = E v h v ~ n ^ Vw'PvPv', (16) 

v vv' 

H' = - v uu '(n u - p u )(h v > - p v >), (17) 

vv' 



Here, we minimize H' by taking p's in the following way 

n = /« \ = ir Vol e = (i9) 

Hu X Vl Xr e -{Ho-liN)/k B T e {E v -y)/k B T _l V 1 
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Then, the thermodynamic potential f2 = —pV can be given as 
n = -k B T\nTr e -("o-/^)/ fefl T 

- -^^'K)K') + ^T^ln(l- e -^'-^ T ). (20) 

vv 1 v' 

We note that Eq. (19) is also determined by the condition 

(or^-J = -V*w(rv> + YV,'-- — ^— — - = 0.(21) 

This condition means that the thermodynamic potential Q is an extremum 
with respect to any change of T, V and \i in a state of thermal equilibrium. 



2.2 HFB Theory 

Now, we present the HFB theory. This theory is reliable for the description of 
the low temperature behavior of the Bose gas. In this theory, the operators 
a and a\ are replaced by the c-number a , a\ w \/T%. Then, the HFB 
Hamiltonian K = H — jiN is given by 

K = (Too-^iVo + iCX 2 

+ ^ E + (22) 

where n A is the expectation value of the non-condensate density of particles, 
and we eliminate the lower power of Nq and terms which are proportional to 
a u and a\. Now, we define new boson operators by Bogoliubov transformation 

c u = u v a v + v u al, (23) 
4 = u u al + v v a v , (24) 
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where u v and v u satisfy u 2 u — v% = 1. Then, we can write the Hamiltonian in 
terms of the new operator c u 

K = (T o - V)N + \v™N 2 Q + ( £ ™' v " - NoV^u v v vl ) 

i/i/Yo 

+ ^2 { e uu'{u u U u i + v v v v >) - N V^,(u u v u > + v v u v i) } 4<v 
w'/0 - 



+ <J S vv >U v V v i + ^oK°i/°'(^^' + V v Vui) \ (44, + C„<V 



(25) 



where e^/ is given as 



= T^-^, + 2(^o +n ^). 



(26) 



We can eliminate the terms proportional to 4 c t' + C ^<V by imposing the 
condition 



e vv iu v v v i + -N Q V™,(u v u v > + v v v v i) = 0, 



(27) 



or 



e^/tV - N V™,v v ' = E v ,u v >5 vvl , 



(28) 
(29) 



where the eigenvalue E u is given as 

Eu = \]el v - (iVoVT) 2 - 
Then, we obtain the diagonal Hamiltonian 

K = (T m -^N Q + l -V™Nl + Y,EAc»- 



(30) 



(31) 



Here, we note that the first two terms in the right hand side of Eq. (31) 
correspond to the GP equation, and the third term in the right hand side of 
Eq. (25) represents small correction terms, and therefore we can ignore the 
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third term. From Eq. (31), we can also obtain the thermodynamic potential 

n = -k B T\nTre- k/kBT 

* (T 00 - fi)N + ±Otf + k sT Yl ln (! - e- E » /kBT ). (32) 

The chemical potential is given by the condition 







cW ' 

= T 00 -/i + C^o + 2^<V; °, (33) 
where we ignore the anomalous average. 

We note that Eqs. (31) and (33) are similar to results of the HFB approach 
to the Popov approximation ||. In the next section, we will see the ground 
state stability of the interacting Bose gas using the GP equation. 

3 Ground State Stability 

In this section, we study the ground state properties of the Bose gas interact- 
ing through the Van der Waals type — CqT" 6 and the Car -3 interactions. We 
note that we must introduce a cut off at r = 2(R), since these interactions 
diverge at r = 0. Here, we assume the hard core potential inside the atomic 
radius (R). This is chosen to be the exponential or Cur~ 12 potential [12| . 



Now, we study the ground state stability of Bose gas for the condensate state 
at T < T c . In this case, the HFB theory is reliable for the investigation of the 
condensate state. To investigate the ground state stability, we employ a vari- 
ational method. It is a good approximate scheme to study the ground state 
stability. We assume the following Gaussian wave function for the ground 
state 
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where a represents the variational parameter. This choice is natural when 
we take the noninteracting limit. 

First, we consider the Van der Waals interaction. Here, we make comments 
on the scattering length a for this interaction. This can be analytically given 

by m 



T(3/4) /mC 6 V 
2v^r(5/4) V 



1 — tan 



7T 



(35) 



where <fi is the semiclassical phase calculated at zero energy from the classical 
turning point to infinity. This phase depends on the repulsive hard core 
potential. Substituting Eq. (34) into Eq. (31), we can obtain the ground 
state energy E gfj 

3 I 1 < 2 . -2\ -6 



2 Nhuj^-(a 2 + a-')-g 6 a- b j, (36) 
where we rewrite a as a = cr/ciho in units of the harmonic oscillator length 



ciho = ^/h/muj and J 6 denotes the dimensionless integral 

= ,J" dss ^]^ e -^0J^. (37) 



The dimensionless coupling constant g$ is defined by 

16iVC 6 / 6 
® Q "&ixhuja\ 



(3* 
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Table 1. We plot the integral values 
of I 3 and Iq for several cases of 2(R) 



2{R) 




/s 


h 


1.0 x 10" 


-4 


1.295 


1.306xl0 10 


2.0 x 10- 


-4 


1.187 


1.632xl0 9 


3.0 x 10- 


-4 


1.123 


4.835x10 s 


4.0 x 10- 


-4 


1.078 


2.040 x10 s 


5.0 x 10- 


-4 


1.043 


1.044x10 s 


6.0 x 10- 


-4 


1.015 


6.044xl0 7 


7.0 x 10- 


-4 


0.991 


3.806xl0 7 


8.0 x 10- 


-4 


0.970 


2.550xl0 7 


9.0 x 10" 


-4 


0.951 


1.791 xlO 7 


1.0 x 10- 


-3 


0.935 


1.306xl0 7 



In Table 1, we plot the values of Iq as a function of the dimensionless param- 
eter 2(R) = 2{R)/aho- We note that 2(R) and a^ are of the order of A and 
/im, respectively. Therefore, 2(R) is of the order of 10~ 4 . 
Now, we look for the critical coupling constant g§. This is given by the 
inflexion point of E g6 

dE d 2 E 
da da 2 

at a = a c and = g^. From Eq. (36), we obtain 

96 = ^- 0.0417. (40) 



Next, we consider the C^r 3 interaction [18"]. We can write down the ground 
state energy E gz 

E 93 = -Nhuii* 2 + a- 2 ) + -5-|/ 3 , 

= ^Nhu |i(a 2 + a- 2 ) + ^a- 3 }, (41) 

where /3 is also a dimensionless integral defined as 

poo pco 1 

h = 2 1 dss2e ~"L m dtt ' e ~"w^m (42) 



n 



In Table 1, we also plot the values of I3 as a function of 2(R). The dimen 
sionless coupling constant #3 is defined by 

16iVC 3 / 3 



93 



3irhua 3 ho 



(43) 



In this case, the ground state energy has the same property as the contact 



interaction p9| . We can obtain the critical coupling constant #3 ~ —0.178. 
Finally, we consider the C%r~ 3 — C?,r~ 6 interaction. Since the Van der Waals 
interaction is weak, we always assume |g 3 | > g§ > 0. The ground state energy 
Eg is given by 



E n 



^Nhcu l^(a 2 + a' 2 ) + g 3 a- 3 



-6 



(44) 



In Fig. 1, we show the ground state energy in units of3/2Nhu as a function 




Figure 1: We show the ground state energy Eq. (44) in units of 3/2Ntkv as a 
function of a for several values of the parameter g = gs/g^. The solid line is 
drawn with g = g c . The dashed and dot-dashed lines are g > g c and g < g c , 
respectively. 



of a for several values of the parameter g = g$/ ' g§ which is given as 

C3I3 a 



9 = T^fat. (45) 
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As can be seen, for g > g c , the ground state energy always has a minimum. 
Therefore, the ground state is stable. On the other hand, for g < g c , a 
minimum disappears. In this case, the ground state is unstable, therefore, 
BEC collapses. We can also calculate the critical value g c = gzjg%. It is given 
by 



= 8(7° 



18^6 



4a c , 

4 



(46) 
(47) 



Coupling constants g^ and g e are related to each other through Eqs. (46) 
and (47). We note that the parameter a c must be a c > 0.66874, since we 
assume g 6 > 0. Eliminating a c from Eqs. (46) and (47), we obtain 

5/4 4 fl + y/1 + 360^ 1/4 



9s 



8 + + 360^ 6 

9 V 10 



9 



10 



(48) 



In Fig. 2, we show the phase diagram of coupling constants #3 and ge- The 



93 



0.5 
0.4 
0.3 
0.2 
0.1 


-0.1 
-0.2 



stable 






unstable 



0.05 



0.1 



0.15 



0.2 



9s 



Figure 2: We show the phase diagram of coupling constants g^ and g$. The 
curve corresponds to the critical curve given by Eq. (48). The parameter 
runs 0.66874 < a c < 1 in this diagram. Above the critical curve, the ground 
state energy is stable, and below the critical curve, it is unstable. 



curve corresponds to the critical curve given by Eq. (48). Above the critical 
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curve, the ground state energy is stable, and below the critical curve, it is 
unstable. Therefore, it is seen that the BEC collapses below the critical curve 
in the phase diagram. 

Now, we estimate the coupling constant g for alkali atoms. Coefficients C3 
and C*6 are numerically given by Ref. |20| for alkali atoms. 



Table 2. We plot the ratio \C 3 \/C 6 x 10 3 between S and 
Pa andP7r state for alkali atoms in units of the atomic 
units (a.u.). 





Li 


Na 


K 


Rb 


Cs 


S-Pa 
S-Pn 


5.329 
3.298 


2.995 
2.325 


1.845 
1.377 


1.527 
1.144 


1.205 
0.885 



In Table 2, we plot the ratio \C^\/Cq for alkali atoms in units of the atomic 
units (a.u.). As can be seen from Table 2, \C^\/Cq is of the order of 10~ 3 
a.u.. The ratio I^/Iq and the harmonic oscillator length a^o are of the order 
of 10~ 7 ~ 10~ 8 and 10 4 a.u., respectively. Thus, the coupling constant g is 
of the order of 10 ~ 10 2 . In Fig. 3, we show the critical coupling constant 

2 

1.5 
1 

0.5 


-0.5 
-1 

0.8 1 1.2 1.4 1.6 1.8 2 

Oc 

Figure 3: We show the critical coupling constant g c as a function of <r c . 
Above the curve, the ground state energy is stable, and below the curve, it 
is unstable. 

g c as a function of a c . Above the curve, the ground state is stable, and be- 
low the curve, it is unstable. As can be seen from Fig. 3, for g > 0, the 
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ground state of alkali atoms are always stable in our calculation. On the 
other hand, for g < 0, we must carefully consider the ground state stability, 
since it depends on o c . We can show that the ground state is stable when 
cr c < 0.756 for g = —10, and when a c < 0.690 for g = — 10 2 . For real 
atoms, the value of a c is expected to be the same as the 5-function type in- 
teraction a c = 0.669 [fl9[| . Therefore, the ground state is also stable for g < 0. 

4 Conclusions 

We have studied the ground state stability of interacting Bose-Einstein gas 
with an arbitrary power low interaction. In particular, we have considered 
the Van der Waals — Cqt~ 6 plus C^r -3 interactions. Using the Gaussian vari- 
ational function, we have obtained the ground state energy. It is shown that 
the ground state stability depends on the ratio of the coupling constants g^ 
and g$, then the critical coupling constant is obtained by the inflexion point 
of the ground state energy. For g > g c , we can always produce stable BEC. 
On the other hand, for g < g c , BEC collapses. Here, the ground state sta- 
bility mainly depends on the C^r~ 3 term, since atoms become sensitive to 
this interaction when the density of atoms becomes high at low temperature. 
Therefore, we can understand that the role of this interaction is important 
for the stability of the ground state at low temperature. 
Next, we have obtained the phase diagram of coupling constants g% and ge- 
From the phase diagram, we have shown that there exists a stable region for 
Bose-Einstein gas. It is useful because we can classify atoms which become 
BEC. 

Finally, we have argued the validity of these results by using the numerical 
values of C3 and Cq for alkali atoms. We have shown that the ground state 
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of alkali atoms are always stable. 

Our results show that we can reliably understand Bose-Einstein gas interact- 
ing through the realistic interaction beyond the GP theory. 
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